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ABSTRACT 

Mass segregation is observed in many star clusters, including several that are less 
than a few Myr old. Timescale arguments are frequently used to argue that these 
clusters must be displaying primordial segregation, because they are too young to be 
dynamically relaxed. Looking at this argument from the other side, the youth of these 
clusters and the limited time available to mix spatially distinct populations of stars 
can provide constraints on the amount of initial segregation that is consistent with 
current observations. We present n-body experiments testing this idea, and discuss 
the implications of our results for theories of star formation. For system ages less 
than a few crossing times, we show that star formation scenarios predicting general 
primordial mass segregation are inconsistent with observed segregation levels. 
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1 INTRODUCTION 

Many star clusters show some degree of radial mass segrega- 
tion, with massive stars over-represented in the inner region 
of the cluster. This is expected in bound clusters that are 
older than a dynamical relaxation time, such as the Arches 
near the galactic center, which shows an almost-flat mass 
function near the center ( |Kim et aL 2006). Theoretical expla- 
nations and numerical explorations of this segregation pro- 
cess are numerous (e.g|Spitzer|p69l|Spitzer fc Shull|[l975 



Khalisi et al.||2007| [Portegies Zwart et al.|2007[ 

Younger clusters, with ages less than a few Myr, are not 
expected to display the same degree of segregation, as they 
are only a few crossing times old and should not have had 
time to segregate dynamically. Despite this, there are several 
examples of young (or even embedded) clusters exhibiting 
mass segregation. Examples with ages less than ~ 5 Myr in- 
clude: Mon R2 ( |Carpenter et al.|1997|>; IC 1805 ( |Sagar et al" 
1988]); NGC 1893 ||Sharma et al.||2007l>; NGC 6530 (McNa- 



mara fc Sekiguchi|1986| ); NGC 6231 ( |Raboud fc Mermilliod' 
1998[); and the Or ion Nebula Cluster (ONC) ( |Hillenbrand 
fc Hartmann|1998| . Some of these arguably show a general 
trend of segregation, with more massive stars progressively 
more centrally located. Others, such as Mon R2, show little 
evidence or this type of segregation, yet still have a central 
concentration of the most massive stars. The youth of these 
clusters relative to their relaxation time makes the process 
of dynamical segregation less likely, and this timescale argu- 
ment has been used as evidence that primordial segregation 



has played a role (Hillenbrand & Hartmann 1998 Bonnell 
fc Davies|1998l|Raboud fc Mermilliod|1998| ). 



Most numerical explorations of dynamical mass segre- 
gation have focussed on long timescales (lOs-lOOs of cross- 
ing times, e.g. |Khalisi et al.|2007| . |Bonnell fc Davies| (j!998 \ 



explored the short-timescale n-body evolution of a variety 
of initial cluster setups, concluding that the formation of 
a Trapeziumesque system in the center of an ONC-scale 
cluster is unlikely to be realized by two-body relaxation 
alone, and that the most massive stars in Orion would have 
to have formed within the inner ~ 25% of the cluster to 
find themselves in the center at the present day. [McMil- 
lan et al. (12007) performed n-body experiments of merging 



sub-clusters, showing that the merger products can exhibit 
a more advanced state of segregation than their age would 
suggest, somewhat ameliorating the timescale problem as 
the memory of rapid dynamical evolution in small- n clus- 
ters is apparently retained when they merge. Nonetheless, 
it remains unclear how important primordial segregation is 
relative to later dynamical evolution. 

While mass segregation is a robust and well studied re- 
sult of dynamical interactions in a cluster, the usual starting 
state for these explorations is an initially unsegregated clus- 
ter. It is unclear how appropriate such initial conditions are 
for very young clusters. Theories of massive star formation 
(e.g. |Bonnell et"aT1|2001| |McKee fc fanl|2003[ ) can easily 
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be interpreted to imply that the mass of a star depends 
on its location in such a way that the structure of a cluster 
when gravitational dynamics become dominant is highly seg- 
regated. In this case, mass segregation is not the eventual 
result of a cluster trying to achieve equipartition, but the 
(perhaps unstable) starting point. In this paper we approach 
mass segregation from the other side, exploring the evolu- 
tion of very young clusters from a primordially segregated 
state. We present n-body experiments covering a range of 
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initial cluster setups and compare the distributions of stars 
to observations of clusters only a few crossing times old, and 
we discuss the implications of our results for theories of star 
formation. 



2 NUMERICAL DETAILS 

To explore this problem, we performed a suite of TV-body 
experiments using the code NBODY6 (|Aarseth||2003|). The 



code uses a fourth-order Hermite integrator to advance the 
self-gravitating system. The energy error AE/Eo over the 
course of each simulation was held under 10 -4 . 

The stars in the simulations have masses in the range 
0.08 - 50.0 MfD drawn from the IMF given by Kroupa 



(2001) 
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0.08 < m/M Q < 0.5, 
0.5 < m/M Q < 50.0 



(1) 



Throughout, we neglect the effects of an initial binary 
population or the remnants of the cluster's gas reservoir. In 
order to begin with a mass-segregated cluster, we set up the 
positions of the stars in a mass-dependent fashion. Each star 
is assigned a uniformly distributed location within a sphere 
of radius 



(2) 



where M is the mass of the star being initialised and M m i„ 
is the minimum mass in the cluster, the index q controls the 
level of segregation; with q = all stars have the same max- 
imum radius independent of mass, while q = 1 yields a very 
strict radial ordering by preventing massive stars from being 
located in the cluster outskirts. At this point the stars are 
spatially arranged with a radial mass (or number) density 
that does not correspond to any standard scenario. Working 
from r = outward, each star in turn has its radius adjusted 
to correspond to some underlying distribution. The first step 
thus effectively ranks the stars radially as a function of their 
mass, while the second step adjusts the initial conditions to 
a physically motivated radial profile. 

Sorting the stars radially by mass has a drastic effect on 
the dynamical structure of the cluster. The mass profile of a 
heavily mass segregated cluster with an isothermal number 
density profile is much more centrally condensed than an 
unsegregated cluster. To encompass a wide range of initial 
parameters, we perform two sets of simulations. In the first 
set, we adjust the radii of the bodies so that the number den- 
sity is uniform with radius, n(r) = n const . In the second set, 
we adjust the radii so that the mass density is an isother- 
mal sphere, p(r) oc r~ 2 . This segregatio n scheme pr e serves 



the overall cluster IMF. We note that Subr et al 



(2008) 



recently devised a method for generating stable, virialized 
clusters with a parameterized degree of mass segregation, 
and we defer consideration of this method to the discussion 
section. 

In each set, we performed three series of 100 experi- 
ments, each with a different level of initial mass segregation: 

Series A - no initial segregation: we set q — in equation 
[2] so that the radial distribution is independent of mass. 

Series B - moderate initial segregation: we set q — 0.15 
in equation [2] At this value, stars more massive than ~ 




Figure 1. Evolution of the half-mass radius r/, m for each of the 
uniform number density series, relative to the initial half-mass ra- 
dius ro. The median value of all the simulations is plotted, with 
error bars showing the upper and lower 95th percentile values. 
Horizontal offsets separate the error bars. The different lines cor- 
respond to the following subsets of the system: thick: all stars; 
solid: 0.08 ^ m/A/ < 1.0; dashed: 1.0 ^ m/M Q < 5.0; dotted- 
5.0 ^ m/M Q 50.0. 



5 Mq are restricted to lie inside half the cluster radius. 
After adjustment to an underlying distribution, this changes 
somewhat. 

Series C - heavy initial segregation: we set q = 1 in equa- 
tion [2] This value of n effectively sorts the stars almost 
strictly by mass, with the most massive stars in the clus- 
ter center. 

Each cluster begins with 1000 stars, and we run each 
cluster for 6 crossing times. We scale all of the clusters so 
that the initial half-mass radius rum = 0.5 pc, the rms ve- 
locity in the clusters is a ~ 1.5 km s _1 for each series, and 
the initial crossing time r c = 2rh m /a « 0.65 Myr. With 
this scaling we are simulating the clusters for more than 3.5 
Myr, encompassing a more-than-sufficient time to compare 
to young clusters such as the ONC. 



3 UNIFORM NUMBER DENSITY RESULTS 

This set of simulations began from an initially uniform num- 
ber density, upon which mass segregation was imposed with 
no concern for what this does to the mass distribution. All 
of the simulations in this set begin in virial equilibrium. In 
our analysis of the simulations, we split the stars into three 
groups: low-mass stars, 0.08 < M*/Mq < 1.0; moderate- 
mass stars, 1.0 ^ M*/Mq < 5.0; and high-mass stars, 
5.0 < M„/M & < 50.0. 



3.1 Radial Evolution 

We begin by considering the overall evolution of the clusters 
in this set, measured by the half-mass radii of the different 
mass components. In Figure [I] we show the evolution of the 
half-mass radius for each series, broken into mass bins and 
over the whole cluster. Plotted is the median value over all 
the simulations in a series. Examining the trends in these 
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Figure 2. Demonstration of our reduction scheme, for stars more 
massive than 5 Mq in series A at t = 6t c . The cumulative number 
distribution with radius for each of 100 runs are shown as black 
lines. The dark shaded region shows the 95% confidence interval of 
the location of the underlying distribution's median, found via a 
bootstrap resampling of the simulation data. The outer lines mark 
the uper and lower 95th percentiles of the simulated distribution, 
to give an idea of the scatter in the ensemble. Further figures in 
this work will show only these lines and the median. 



quantities for series Al (initially unsegregated) , it is appar- 
ent that after only a few crossing times some degree of mass 
segregation is setting in. While the overall half-mass radius 
remains fairly constant, massive stars have moved inwards 
so that after 6 crossing times their r^ m is ~30% of its ini- 
tial value. This evolution occurs as the least-massive stars, 
and to a lesser extent the moderate-mass stars, are driven 
outward. This degree of rapid segregation is con sistent with 
Bonnell fc Daviesf |l998| and |Subr eiH] ( |2008 \ . 

Comparing the evolution of the half-mass radii for se- 
ries Bl (moderate initial segregation) with series Al, the 
trends are similar. With the exception of the high-mass bin, 
the changes in rum for each component in series Bl closely 
match the unsegregated case but offset by the initial segrega- 
tion. The inward evolution of r^m for massive stars appears 
to saturate at ~ 30% of the total cluster value, so that after 
t = 5 t c series Al and Bl have similar values. The overall 
half-mass radius is nearly constant over the course of the 
simulations. 

The overall cluster evolution of series CI (heavy ini- 
tial segregation) is qualitatively different than the other two 
uniform number density experiments, as the cluster expands 
from its initial configuration. Rather than the massive stars 
moving inwards as in series Al and Bl, the massive popula- 
tion remains radially stationary. The low-mass stars adjust 
their location quickly outwards between 1 and 2 crossing 
times, and the moderate-mass stars follow a general out- 
ward trend with time. 



3.2 Cumulative Distributions 

For each mass range, we consider the cumulative number 
distribution of stars. Comparing these distributions gives a 
measure of the mass segregation in the systems. We treat the 
results of each series of initial conditions as an ensemble, and 
attempt to recover the underlying distribution for each mass 
bin at different times. 



In Figure [2] we illustrate this process. In this example, 
from series A at t = 6 r c , we plot for all 100 experiments 
the cumulative fraction of high-mass stars as solid lines. This 
population of stars displays the largest spread in the the fi- 
nal distribution within one series, due to the relatively low- 
number statistics (the mean number of stars in this mass 
range is 11, with standard deviation ~ 3). The spread of 
the cumulative distribution function is asymmetrical, so we 
use the median as the characteristic value for the ensem- 
ble of simulations. The dark shaded region shows the 95% 
confidence interval of the underlying distribution's median, 
obtained by a bootstrap resampling of the simulation data. 
The thick lines mark the 95th percentiles of the simulated 
ensemble, and give an indication of the spread in simulated 
values. In the rest of this discussion we plot only these quan- 
tities, rather than the results of each individual run. 



3.2.1 Series A 1- Uniform number density, no segregation 

Figure [3] shows the time evolution of the distributions of the 
three mass components for all three series in the uniform set, 
with the left column showing series Al. In the four panels we 
show the cumulative distribution of each mass component at 
times t = 0, 1,3 and 6 r c . The radius in each plot is shown 
as a fraction of the half-mass radius at that time, rather 
than at the initial time. We make this choice to show what 
an observer might see looking at a cluster. With knowledge 
only of the current state, measuring radius relative to an 
unknown and irrelevant initial mass distribution would be 
impossible. 

At the beginning of the cluster evolution, the stars 
in each mass bin are uniformly distributed throughout the 
cluster, evidenced by the identical form of their initial cu- 
mulative distribution functions. After 1 crossing time the 
shape of the cluster has changed somewhat, and even at 
this early time there is a hint of mass segregation of the 
high-mass stars. The distributions of the different compo- 
nents are mostly still overlapping, however. By 3 crossing 
times the segregation of the massive stars has become more 
pronounced, while the moderate- and low-mass stars show a 
small separation in their distribution but remain quite sim- 
ilar. Between 3 and 6 crossing times the separation between 
the low- and moderate-mass distributions becomes clearer, 
though their medians are nearly coincident with the other's 
lower and upper 95th percentile respectively. From t — 3r c 
onward, the high-mass distribution is completely separated 
with a high degree of confidence from the other two com- 
ponents in the inner regions of the cluster, and the median 
value is separated from the extremes of their distributions 
at all radii. 



3.2.2 Series Bl-Uniform number density, moderate 
segregation 

The middle column of Figure [3] shows the cumulative distri- 
bution plots for series Bl. In this model the massive stars are 
all initially placed within the half-mass radius, but there is 
still a radial overlap among radial distributions of the three 
mass components (i.e. low-mass stars are present in the in- 
nermost region of the cluster; this is in contrast to series CI). 
Because of the central concentration of the massive stars, 
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Figure 3. Evolution of the cumulative number distribution of three mass components for the series with uniform number distributions. 
Columns are labeled with the series they represent, and rows are labeled with time increasing downwards. Filled region: 95% confidence 
limit on the median values of the distribution; lines: 95th percentile limits of the distribution. 



the half-mass radius is a smaller fraction of the total cluster 
radius. The evolution is very similar to that seen in series 
Al, with the separation between the low- and high-mass 
distributions getting clearer with increasing time, but with 
more advanced levels of segregation throughout. The main 
difference between this series and the unsegregated case is 
the separation between the low- and moderate-mass stars. 
The initial difference in these distributions is not erased over 
the course of the simulations. This is best seen by compar- 
ing the results from series Al and Bl at t = 3r c . In the 
unsegregated runs, the medians of the low- and moderate- 



mass distributions are marginally separate, but the lower 
95th percentile of each distribution is nearly the same. In 
contrast, the lower-mass components are still clearly sepa- 
rated at this time in series Bl, with the respective upper and 
lower 95th percentiles nearly coincident for much of the ra- 
dial range. The difference between these runs is less distinct 
after 6 crossing times. 
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3.2.3 Series Cl-Uniform number density, heavy 
segregation 

The cumulative distributions for this series are shown in 
the right-hand column of Figure [3] The extreme nature of 
the initial segregation is apparent in the lack of radial over- 
lap among the three populations; the moderate-mass stars 
only begin to appear outside of the high-mass stars, and the 
low-mass stars are almost exclusively located outside the 
half-mass radius. Up through t = 3r c the evolution of the 
distributions is mostly dominated by the outward movement 
of the low-mass stars, with the separation between the medi- 
ans of the different distributions remaining clear at all radii, 
and the extremes of the distributions distinguished inwards 
of rhm- By t — 6t c the moderate- and high-mass distribu- 
tions have moved closer together, though inwards of ~ 2n lm 
they remain clearly separated. 

The point of this series is to test the efficiency of mixing 
from such an extreme setup. If this case showed thorough 
mixing after only a few crossing times, than it would seem 
that observations of mass segregation in young clusters hold 
little value in telling us where stars of different masses form. 
However, the clear signal of segregation is retained through 
at least 6 crossing times, corresponding to several Myr for a 
parsec-scale cluster. 
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Figure 4. Evolution of the half-mass radius r^m f° r each of the 
isothermal mass density series, relative to the initial half-mass ra- 
dius ro. The median value of all the simulations is plotted, with 
error bars showing the upper and lower 95th percentile values. 
Horizontal offsets separate the error bars. The different lines cor- 
respond to the following subsets of the system: thick— all stars; 
solid: 0.08 < tti/Mq < 1.0; dashed: 1.0 < m/M Q < 5.0; dotted: 
5.0 ^ m/M Q ^ 50.0. 



4 ISOTHERMAL MASS DENSITY RESULTS 

In these experiments the radii of the stars are adjusted so 
that the underlying mass density is an isothermal sphere, 
p(r) oc r~ 2 . Thus the mass profiles and total radial extent 
of each of the three series are identical, with the outer cluster 
radius at 2rh m , though the number densities are much dif- 
ferent. The non-segregated case has a number density that 
is also isothermal, while the segregated cases have shallower 
number density profiles that are radially-dependent . This 
set of experiments can be thought of as representative of 
clusters forming from the same progenitor molecular cloud, 
with different star formation mechanisms. The results pre- 
sented all began in a state of virial equilibrium, though we 
also ran these as slightly sub- virial systems with a virial pa- 
rameter q — Ekin/\E po t\ — 0.4. Differences are noted where 
appropriate below. 



potential energy increases relative to an unsegregated clus- 
ter, and the velocity dispersion consequently increases. The 
tendency for fast-moving stars at the cluster outskirts to 
drift away from the cluster is exacerbated somewhat by not 
imposing a radial dependence on the velocity dispersion, as 
in a Plummer or King profile. Though this would bolster the 
stability of the radial mass profile, observations tend to favor 
single-dispersion descriptions with little evidence for radial 



dependence in young clusters (for instance in the ONC, Hil- 
lenbrand fc Hartmann|1998 1. 



In series C2 (an as in series B2 and CI) there is a rapid 
outward adjustment of the low-mass stars between 1 and 
2 t c , after which the evolution of this series is very similar 
to the uniform number density case. A moderate increase 
of the overall cluster radius takes place as the massive star 
distribution remains largely stationary. 



4.1 Radial Evolution 

We again begin by examining the half-mass radii for this 
set, shown in Figure [4] The top panel shows series A2. This 
series is quite similar to the uniform number density case, 
though with a small degree of overall cluster expansion over 
6 crossing times. The same trend of inward migration for 
massive stars with outward movement of the low mass stars 
is seen, and to a similar extent as in series Al. 

The main difference between the half-mass radii evolu- 
tion in series B2 versus Bl is the early outward movement 
of the low mass stars between 1 and 2 r c . The root cause of 
this, which is also seen in series C2 and was present to some 
extent in series CI, is the assumption of a mass-independent, 
isothermal velocity distribution. When mass is concentrated 
to the center of the cluster, the magnitude of the system's 



4.2 Cumulative Distributions 

In Figure [5] we show the evolution of the mass distribution 
for this set in the same fashion as for the uniform number 
density results. 



4-2.1 Series A2~Isothermal mass density, no segregation 

The isothermal mass distribution for this set is reflected in 
the linear growth of the cumulative fraction for each mass 
component. After 1 crossing time there is no real separation 
between the components. By 3 crossing times the high-mass 
distribution is beginning to separate from the lower mass 
stars, which continues through t = 6r c . As in the previous 
setup with no segregation, only marginal separation between 
the distributions of the low- and moderate-mass stars is seen. 
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Figure 5. Evolution of the cumulative number distribution of three mass components for the series with isothermal mass distributions. 
Columns are labeled with the series they represent, and rows are labeled with time increasing downwards. Filled region: 95% confidence 
limit on the median values of the distribution; lines: upper and lower 95th percentiles of the distribution. 



4-2.2 Series B2-Isothermal mass density, moderate 
segregation 



The middle column of Figure [5] shows the cumulative distri- 
bution evolution for this series. As in the uniform number 
density case, the massive stars are initially all within the 
half-mass radius, but there is still radial overlap among the 
components. The initial separation between the components 
never gets erased through 6r c , though by the end of the runs 
the extremes of the distributions are overlapping. The me- 



dians, however, remain clearly separated at all radii through 
3 crossing times, and within 2rh m at 6 crossing times. 

When evolved from slightly sub-virial initial conditions, 
the evolution of this series is only mildly affected. The period 
between t = and t = lr c is marked by overall contraction 
of the cluster as it seeks equilibrium. This contraction is seen 
in all three mass components. From that point onwards, the 
behavior of the sub-virial system is as described for the virial 
case, with expansion of the low-mass component and slow 
growth of the overall half-mass radius. The changes to the 
cumulative distribution plots are minimal. 
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Figure 6. Evolution of the cumulative number distribution of 
three mass components for series D, with an isothermal mass dis- 
tribution and initial segregation of only the massive stars. Times 
are labeled on the figure. Filled region: 95% confidence limit on 
the median values of the distribution; lines: upper and lower 95th 
percentiles of the distribution. 



4-2.3 Series C2-Isothermal mass density, heavy 
segregation 

Like the other segregation levels, the evolution of this se- 
ries shows a close similarity to its uniform number density 
counterpart. The initial conditions again show an almost- 
total amount of segregation, with the moderate-mass stars 
only appearing outside of the massive stars, and the low- 
mass stars excluded from inside the half-mass radius. The 
separation between the moderate- and high-mass distribu- 
tions grows less distinct with time and increasing radius, but 
they remain separate at a level comparable to series B2. The 
low-mass distribution is separate from the higher mass stars 
even at the 95th percentile level out to ~ 2r^ m at all times. 

As with the moderate segregation case, running this 
system sub-virially in an effort to stem the expansion of low- 
mass stars has a minimal impact on the overall evolution, 
with the only real difference during the first crossing time 
as the cluster contracts to ~ 70% of its initial size before 
evolving as the virial case does. The cumulative distribution 
plots are likewise almost unaffected. 



4-2.4 Series D-Isothermal mass density, only massive 
stars segregated 

If massive stars form via a fundamentally different process 
than low-mass stars (for instance, if stars earlier than type 
B require different conditions to accrete past the obsta- 
cle of their own radiation), it could be the case that only 
these stars would be initially segregated. This idea is sup- 



ported by dynamical studies ( Bonnell fc Davies||1998 1, per- 
haps some observations (e.g. the Trapezium Hillenbrand & 
Hartmarm]|1998| , and is implicit in some theories of mas- 
sive star formation (e.g. the fiducial model of McKee & Tan 
2003). When comparing these results to observations, it will 
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Figure 7. The ratio of the number of low-mass to moderate-mass 
stars within the half-mass radius, Ni ow /N mel i- The median value 
for each series is plotted, with error bars showing the upper and 
lower 95th percentile values. The small horizontal offsets serve to 
separate the error bars, all times are measured at integer crossing 
times. 



mass density scheme to test this scenario, and the results 
are shown in Figure [6] All of the massive stars are initially 
inside the half-mass radius, while the distributions of low- 
and moderate-mass stars are identical and extend through- 
out the cluster. Throughout the run the massive stars remain 
highly segregated. The low- and moderate-mass stars segre- 
gate at a level comparable to the unsegregated cases Al and 
A2. A scenario in which the most massive stars form in a 
compact configuration at a cluster center, as the Trapezium 
appears now and as W3 IRS 5 appears at a younger age 



therefore prove useful to consider a setup in which only the 
massive stars are initially segregated. We use the isothermal 



| |Megeath et al. 2005 I , is difficult to distinguish from ini- 
tially unsegregated initial conditions if one looks only at the 
relative distributions of the low- and moderate-mass stars. 



5 DISCUSSION OF NUMERICAL RESULTS 

Examining the distributions in Figures [3] and [5j the clearest 
discriminant between the different levels of segregation ap- 
pears to be the distribution of the low- and moderate-mass 
stars. The unsegregated cases show only a very marginal 
separation between these distributions, while the runs with 
initial segregation maintain the initial distinction through- 
out the simulations. To put a number to the visual difference 
between these populations, we use the ratio of low-mass to 
moderate-mass stars, Ni ow / N me d, for all stars located inside 
the half-mass radius. While it could be tempting to compare 
these results to observed clusters, this ratio is extremely sen- 
sitive to survey completeness, and the cumulative distribu- 
tion is probably a preferable means of comparison to reality, 
which we examine further below. 

Shown in Figure [7] is the median value of this ratio for 
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each series, with the error bars showing the upper and lower 
95th percentiles. While there is a fairly large spread in this 
quantity, the median values remain distinct at all times. For 
the uniform number density set, the moderate-segregation 
case (series Bl) begins to lose its distinction from the unseg- 
regated case (series Al), also noted in the discussion of the 
cumulative distribution. The isothermal mass density series 
retain their difference for longer times, also evident in the 
cumulative distribution plots. 

A comparison of the cumulative distribution functions 
from the simulations starting from a uniform number density 
(Figure [3} and from an isothermal density distribution (Fig- 
ure]^) shows a striking similarity between the two cases. The 
only clear difference is in the shape of the initially unsegre- 
gated series, with the uniform number density case showing 
more final segregation than the isothermal mass density set. 
This similarity, despite very different initial setups, argues 
for the robustness of the results to variations in the initial 
conditions. 

A source of potential concern is in our parametriza- 
tion of mass segregation. Our scheme sorts masses radially 
according to mass and then stretches the radii to match 
a chosen distribution, with no concern for the underlying 
stability or resultant dynamical structure. One could argue 
that the star formation process itself may be likewise uncar- 
ing when it comes to the initial setup of a cluster; a young 
cluster is not expected to be in a relaxed state immediately 
following the expulsion of gas and its entrance into a stellar- 
dynamically dominated phase. Nonetheless, a comparison 
to an alternative mass segregation scheme would seem war- 
ranted to ensure that our results are not dominantly a re- 
flection of instabilities in our initial conditions. 



Subr et al 



(2008) (hereafter S08) developed a method 



for constructing initial conditions corresponding to dynam- 
ically evolved, segregated clusters. We chose not to adopt 
their method for this work because we are concerned with 
very young clusters and primordial versus dynamical segre- 
gation. Their system sets up the cluster in a stable, virialised 
state, and is parametrized by an index S. With S — they 
have an unsegregated Plummer sphere, and S = 0.5 yields 
a heavily segregated cluster with a central density that is 
approximately isothermal. Using the code PLUMIXQ we 
generated initial conditions corresponding to S — 0.25 and 
S — 0.5 in S08. We then applied the same cumulative dis- 
tribution analysis to these initial conditions as we used with 
our simulations. These distributions are shown in Figure [8] 
Their S = 0.5 has a near-isothermal mass distribution, 
so comparison to our series C2 would be the most equivalent. 
The distributions from S08 bear little relation to our initial 
conditions, which is not unexpected due to the very different 
parametrization than used here. However, after one crossing 
time the distributions appear remarkably similar. Our initial 
conditions, constructed without regard to stability, seem to 
relax within one crossing time to the stable solutions of S08. 
Further evolution confirms this, as seen in a comparison of 
the half-mass radii for series C2 in our Figure [4] with figure 
5 of S08. We therefore see that our initial conditions - with 
very different initial density profiles showing the same trends 



in the evolution of mass segregation - quickly (within one 
crossing time) approach a form corresponding to a dynami- 
cally mature state of segregation. This level of segregation is 
maintained for several crossing times, as shown in our simu- 
lations and those in S08. We conclude that strong primordial 
mass segregation, largely independent of its paramateriza- 
tion, is unlikely to be dynamically erased over the first few 
Myr of cluster evolution. 



6 MODEL LIMITATIONS AND 
SIMPLIFICATIONS 

In order to keep this study compact and tractable, we have 
made several simplifying assumptions. A realistic star clus- 
ter will display some or all of the following complicating 
factors, which would affect these results to varying degrees. 



6.1 Primordial Binarity 

We have assumed that the initial binary fraction is zero, a 
situation at odds with observation and theory, which show 
that binarity is established to some extent during a stellar 
system's pre-main-sequence evolution (e.g. Mathieu 1994 
Bate 20091. Close encounters between binaries and other 



clusters members can alter the two-body relaxation that 
moves the system toward equipartition (and thus mass seg- 
regation), and could perhaps change the rate of mixing in 
these simulations. However, the short timescale and moder- 
ate stellar densities and velocities in these simulations sug- 
gest that the rate of close encounters will be low, and the 



effect of binaries will be minimal. Adams et al. ( 2006 1 tested 



the non-importance of binarity in similarly sized clusters as 
ours, over a longer time period, and verified this argument. 

6.2 Gas Removal 

We begin our integrations from a gas-free state. The details 
of gas removal will, in reality, affect a cluster to the point 



of determining its survival as a bound entity (Baumgardt & 
Kroupa||2007 1 . The parameter space of gas removal details 
is quite large, and not tightly constrained. In n-body simu- 
lations, it is typically modeled by globally lowering the gas 
mass, effectively removing its potential, over some timescale 
(e.g. I Adams et al.|2006||Baumgardt fc Kroupa|2007||Throop 



& Bally 2008). Generally, the effect of gas removal is to 



1 Authored by L. Subr, and available at http://www.astro.uni- 



bonn.de 



leave the stars moving super-virially, with consequent clus- 
ter expansion. Global expansion of the clusters will have the 
qualitative effect of increasing the amount of mixing at large 
radii, as potentially fast-moving stars from the interior catch 
up with the expanding cluster outskirts, while the inner re- 
gion retains whatever segregation it has as the interaction 
rate between stars drops. This is seen to a mild extent in 
our series C runs, both of which display some level of global 
expansion, and which retain their segregation inwards of 2 
Tfom • 



6.3 The Galactic Environment 

These experiments have taken place in a level of isolation 
that can only be found on a computer. The importance, or 
non-importance, of the galactic gravitational field can be 
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estimated by comparing the cluster size to the tidal radius 
rt , which for a cluster of mass M c on a circular galactic orbit 
of radius d we approximate as the Jacobi limit ( Binney & 
Tremaine|1987l ), 

1/3 



n = d 



f M c 
\3M g (d) 



(3) 



with M g (d) the galactic mass interior to the orbit. Our clus- 
ters all have initial half-mass radii of 0.5 pc, with mass 
Mc ~ 500 Mq. In a galaxy similar to the Milky Way, we 
then have the ratio rnm/rt ^ 0.06 for galactic radii d 6 
kpc. |Baumgar dt fc Kroupa| ( |2007| > showed that for rjj m /r t 
less than this value the cluster evolution proceeds as if it 
were isolated, so the effect of a galactic tide is unlikely to 
affect these results for star formation in the solar neigh- 
borhood. Nearer the galactic center, tidal effects could be 
more important. At d = 2 kpc, for instance, rhm/rt ~ 0.15, 
making the assumption of isolation slightly suspect. Thus 
at later times in these simulations the exterior of a cluster 
near the galactic center may start to be lost to the field, but 
for the short timescales considered here the effect should be 
modest. 



6.4 Cluster dumpiness 

Perhaps the most significant approximation made in this 
study is the spherical symmetry of the clusters. While a 
spherical system offers appealing simplicity and a dramatic 
reduction of parameter space, star formation in nature tends 
to occur in a form that may better be approximated by a 
cylinder, a triaxial spheroid, or more complex filamentary or 
fractal geometries. Breaking spherical symmetry (or intro- 
ducing a bulk motion to parts of the cluster) makes analysis 
of quantities binned by radius problematic, as these may 
be an unsatisfactory average of different regions that have 
evolved to dynamically dissimilar ages. 

The build up of a cluster by merging sub-clusters is seen 
in some large-scale hydrodynamic simulations of cluster for- 
mation (Bonnell et al. 20031. The question of mass segre- 



merger clumps is inherited by the merger product, whether 
this segregation is primordial or dynamically realized on the 
shorter timescales associated with lower-n sub-clusters. This 
yields a large, mass-segregated cluster that appears to be 
several dynamical times older than its 'true' age. The most 
relevant aspect of McMillan et al.'s result to this study is 
that the merging process, akin to violent relaxation, does 
not change the level of pre-existing mass segregation in a 
heirarchical system. 

However, the simulations of McMillan et al. are not di- 
rectly comparable to ours. Firstly, their clusters contained 
over an order of magnitude more stars than the ones pre- 
sented here. More importantly, their initially segregated sub- 
clusters were generated by evolving an unsegregated cluster 
in isolation before merging, a process that treats primordial 
segregation as the stable end point of a dynamical process. 
This work treats primordial segregation as an imposed and 
unstable initial condition, at a higher level of segregation 
than theirs. 

It would be worthwhile (though well outside the scope 
of this work) to explore the question of merging sub-clusters 
with the type of segregation we have used, to see if the vi- 
olent relaxation process mixes the mass components more 
thoroughly than the two-body dynamics considered here. If 
so, this process could ameliorate the inconsistency between 
observations and a strong mass-radius dependence in star 
formation mechanisms that we discuss below. However, if 
the segregation of the merged product remains at the dy- 
namically unattainable levels imposed by our initial con- 
ditions, the nature of our results would be unchanged by 
including sub-cluster mergers. 



gation in such a process has seen preliminary exploration 
n-body exploration by |McMillan et aL| ( |2007| ). Interestingly, 
they show that the level of segregation present in the pre- 



7 COMPARISON TO OBSERVATIONS 

Mass segregation is observed to some degree in many clus- 
ters; while there are clear examples like the potentially ex- 
treme cases or the Arches or 30 Doradus, here we are con- 
cerned with clusters in the earliest stages of their dynami- 
cal evolution. Young clusters (less than a few Myr old) dis- 
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playing some form of mass segregati on includ^] Mon R2 



and N GC 2024 ( |Carpenter et al.|1997||Carpenter fc Hodapp 



20081, which show little evidence for general segregation but 
seems to have a central concentration of the most massive 



sources; IC 1805 ( |Sagar et al.||1988| >; NGC 1893 (|Sharma 
et al.|2007 l, showing segregation ab ove ~ 5 Mq which is ar- 



gued to be primordial; NGC 6530 (McNa mara fc Sekiguchi| 
19861); NGC 623 1 < |Raboud fc Mermilliod||1998[); Stock~8 
(jJose et al.||2008| ); and the ONC ( jHillenbrarid fc Hartmann| 
19981). 



Here we compare observations of two clusters to the 
type of segregation seen in the simulations. The clusters we 
consider are: the ONC, around 1 Myr old and the closest site 
of massive star formation at ~ 400 pc (see the discussion of 



distance in Muench et al. 20081; and NGC 6231, a young 
(3-4 Myr) open cluster 1.8 kpc away ( Raboud et al.|[l997 |. 
These clusters are chosen for their youth and the existence 
of published data sets that permit construction of cumula- 
tive fractional distributions for different mass components. 
In Figure [9] we show the cumulative distributions for these 
sources. It should be mentioned that the radius quoted in 
the observations here is the projected radius, rather than the 
three-dimensional radius shown in our plots. While we have 
chosen to retain the full spatial information in our simula- 
tions, there is very little change to the relative shape of the 
distributions. 



Hillenbrand & Hartmann ( 1998 1 discuss mass segrega- 



tion in the ONC based on data from Hillenbrand (19971, 



which we use to construct cumulative distribution functions 
for the same mass components as in our simulations. The 
ONC shows a clear segregation of stars more massive than 

5 Mq compared to the lower-mass stars, with apparent 
separation of the low- and moderate-mass bins inwards of 
around 1 pc. Two-sample Kolmogorov-Smirnov (KS) tests 
of the distributions confirm that the high-mass distribution 
is probably distinct from the low- and moderate-mass com- 
ponents, with p-values of 0.003 and 0.027 respectively. The 
separation between the moderate- and low-mass distribu- 
tions is statistically unconvincing, with a p-value of 0.087. 

Mass segregation in NGC 6231 is considered by Raboud] 

6 Mermilliod (19981. Data is only available for stars more 
massive than 1 Mq, so we show only the moderate- and 
high-mass distributions for this cluster. The difference be- 
tween the distributions of these two populations is confirmed 
with a KS test giving a p-value of 1.4 x 10 -4 . Raboud et al. 
divide the stars for this data set into several mass bins and 
show evidence for increasing levels of segregation with in- 
creasing mass. We note that if the incompleteness in the 
sample is a function of mass only, and is not radially depen- 
dent, than in principle the underlying distribution of stars as 
traced by the cumulative fraction should remain unchanged. 
If this is the case for NGC 6231, then it is an example of a 
cluster with more general mass segregation at a young age. 

Looking only at the difference between high-mass stars 
and the rest of the stellar population, any of our models 
might be thought to show a similar degree of mass segre- 
gation at times greater than a few crossing times. The best 
discriminant between our models is the separation between 



the low- and moderate-mass distribution. This characteristic 
is not seen in the ONC data with any degree of certainty. 
Series B and C, in both the uniform number density and 
isothermal mass density sets, show a clear distinction be- 
tween the high-mass distribution and the other stars. Series 
A, by contrast, shows very little difference between the low- 
and moderate-mass profiles. Series D, beginning with only 
the massive stars segregated, likewise appears to be consis- 
tent with the available observations. 



8 IMPLICATIONS FOR STAR FORMATION 
THEORY 

The state of star formation theory (especially at the high 
end of the mass function) could not be described as qui- 
escent, with two main theories currently at odds over the 
root mechanism of gathering stellar mass from the parent 
cloud, and how this mass is distributed amongst the stars 
to create the initial mass function. The competitive accre- 
tion theory ( jBonnell et al.|[2001[ ) posits that fragmentation 
occurs at a characteristic mass set by the Jeans instability, 
and accretion onto these protostellar seeds leads naturally 
to a the spread in masses, with massive stars being those 
that were born in advantageous locations. The core collapse 
model ( |McKee fc Tan|2003 1 argues for a direct mapping be- 
tween the mass of a collapsing core and the mass of the star 
that emerges- the mass is assembled prior to collapse, not 
during the early dynamics of an embedded cluster. Massive 
cores are supported against collapse by turbulence and the 
high pressure in massive star forming regions. 

Taking both of these theories at face value and at their 
simplest conceptual limit, one would argue for primordial 
mass segregation. In competitive accretion in a spherical sys- 
tem, the richest gas reservoir is deep in the cluster potential. 
Protostars moving slowly through the cluster center would 
accrete more material than on the outskirts, leading to an 
initially segregated cluster. The core collapse model leads to 
the same situation. The fiducial high-mass core presented in 



McKee & Tan ( 2003 1 supposes that the high mass core is 



located well inside the half-mass radius of the parent clump, 
and it is loosely assumed that the core mass traces the clump 
mass to some degree. If stars acquire their mass from a core 
and don't move far during formation, and the core mass 
is tied to the clump density profile, a degree of continuous 
mass segregation follows naturally. In an idealized spherical 
cluster, both theories would thus seem to argue for a clus- 
ter initially well-segregated, with massive stars preferentially 
forming near the cluster center. 

These results would argue that such a low-level view of 
the theories is unlikely to be reflective of an actual mode of 
star formation. Observations of the youngest clusters, com- 
bined with the level of mixing seen in our simulations, do 
not support a simple spherical cluster with a purely radial 
dependence on the final mass of a star. Understanding star 
formation on a cluster-scale appears to require levels of sub- 
structure not reflected in a spherical theory, and perhaps 
only revealed through cluster-scale simulations. For exam- 
ple, the simulation of a non-spherical GMC collapse in |Bon-| 
nell et al. (2008) displays sub-clustering in the stellar posi- 



2 Though see the cautionary note of |Ascenso et al.] j2008| regard- 
ing evidence for mass segregation. 



tions, an arrangement not considered in this work. 

Despite the lack of segregation across all mass ranges, 
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r/pc 



Figure 9. Cumulative fractional distributions for two young star clusters exhibiting mass segregation. Left: the ONC, with data from 
Hillenbrand ( 1997J and at an assumed distance of 400 pc. Coloured lines correspond to the mass bins found throughout the text. Right: 
NGC 6231, with data from [Raboud &; Mermilliod| ljl998 ) , Data is only available for stars more massive than 1 Mq, and we show only 
the moderate- and high-mass bins. 



it is clear that the most massive stars tend to be centrally 
concentrated at an early age. The Trapezium stars are an 
obvious example, as well as Mon R2, NGC 2024 and NGC 
1893. In Stock 8, cumulative distributions of massive and 
low-mass stars are statistically different if the mass cut is 
4 Mq, while they are indistinguishable if the mass cut is 1 
Mq, suggesting that the most massive stars are segregated 
while the moderate and low-mass stars are more evenly dis- 



tributed (Jose et al. 20081. Raboud & Mermilliod (19981 



argue for a similar state for NGC 6231, though some segre- 
gation of the moderate-mass stars is evident. | Allison et al.| 
1 2009 ) apply a different method, minimum spanning trees, 
to the ONC and find no evidence for segregation below 5 
Mq. A promising route to this arrangement is formation 
in small subclusters that either form segregated or quickly 
dynamically segregate, and merge leaving behind a subclus- 
ter of massive stars in the 'halo' of a less-segregated cluster. 
This process has been seen in cluster-scale simulations and 
has seen preliminary n-body exploration by McMillan et al.| 



1 2007), though more thorough experiments to test the degree 



of post-merger segregation have not yet been performed. 

The alternative to merging subclusters is a fundamen- 
tally different mode of formation for high-mass stars, or spe- 
cial conditions found only at the centers of clusters that lead 
to high-mass formation. This would remove the continuous 
segregation across all masses and leave only the most mas- 
sive stars initially segregated. Observations of young, em- 
bedded massive groupings like W3 IRS 5 (Megeath et al 



2005) might suggest that the tight clustering of massive stars 
is indeed a primordial effect. The results of our series D, 
testing a qualitatively similar scenario, are consistent with 
observations of older (but still young) clusters. 

In summary, dynamical mixing in young clusters is in- 
sufficient to erase the signature of a continuous primordial 
mass segregation that affects stars of all masses. Initially un- 
segregated initial conditions evolve to a state that broadly 
matches the level of segregation observed, but as shown pre- 



viously (Bonnell & Davies 19981 can not account for the 



presence of dense subclusters like the Trapezium. A star- 
formation scenario in which only the most massive stars are 
primordially segregated is consistent with observations, and 



offers a way to account for compact groups of young, massive 
stars. 
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